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ABSTRACT 

We investigate the properties of clusters of galaxies in two cosmological models 
using N-body simulations and the Press-Schecter (PS) theory. In the first model, the 
initial power spectrum of density fluctuations is in the form P{k) oc fc~^ at wavelengths 
A < 120/i~^ Mpc. In the second model, the initial linear power spectrum of density 
fluctuations contains a feature (bump) at wavelengths A 30 — 60ft.~^ Mpc which 
correspond to the scale of superclusters of galaxies. We examine the mass function, 
peculiar velocities, the power spectrum and the correlation function of clusters in both 
models for different values of the density parameter SIq and erg (the rms fluctuation 
on the 8/i~^Mpc scale). The results are compared with observations. We find that in 
many aspects the power spectrum of density fiuctuations in the model (2) fits the 
observed data better than the simple power law model (1). In the first model, the 
mass function and peculiar velocities of clusters are consistent with observations only 
if < 0.6. In the second model, the permitted region in the (JIq, cts) plane is larger. 
In this model, the power spectrum of clusters is in good agreement with the observed 
power spectrum of the APM clusters. This model predicts that there is a bump in 
the correlation function of clusters at separations r ^ 20 — ibh~^ Mpc. In the future, 
accurate measurements of the cluster correlation function at these distances can serve 
as a discriminating test for this model. 

We examine the linear theory predictions for the peculiar velocities of peaks in 
the Gaussian field and compare these to the peculiar velocities of clusters in N-body 
simulations. We determine the clusters as the maxima of the density field smoothed 
on the scale R ^ 1.5/i~^ Mpc and define their peculiar velocities using the same 
smoothing scale as for the density field. The numerical results show that in this case 
the rms peculiar velocities of clusters increase with cluster richness. The rms peculiar 
velocity of small clusters is similar to the linear theory expectations, while the rms 
peculiar velocity of rich clusters is higher than that predicted in the linear theory 
(~ 18% for clusters with a mean intercluster separation dd = ZQh~^ Mpc). 

Key words: cosmology: theory - large-scale structure of Universe, cosmology: theory 
- dark matter, galaxies: clusters 



1 INTRODUCTION 

The formation of structure in our Universe is one of the most 
fascinating problems in cosmology. Usually we believe that 
galaxies and clusters of galaxies have developed by gravita- 
tional instability out of small inhomogeneities of the early 
Universe. The initial field of density fluctuations (5(x, t) can 
be decomposed into its Fourier components S\^{t) and ex- 
pressed in terms of the power spectrum P{k) — {|5kP). 

Figure 1 shows the observed power spectra derived from 
the distribution of galaxies in the APM, Stromlo-APM and 
SSRS2-hCfA2 surveys (Baugh & Efstathiou 1993; Tadros & 
Efstathiou 1996; Costa et al. 1994). The power spectrum 
of the galaxy distribution in the Stromlo-APM redshift sur- 



vey peaks at the wavenumber k — 0.052/i Mpc""*^ (or at the 
wavelength A — VlQh~^ Mpc). A similar peak in the one- 
dimensional power spectrum of a deep pencil-beam survey 
was detected by Broadhurst et al. (1990) and in the two- 
dimensional power spectrum of the Las Campanas redshift 
survey by Landy et al. (1996). Available data, however, are 
insufficient to say whether the peak in the Stromlo-APM 
survey reflects a real feature in the galaxy distribution. It is 
likely that the decline in the power spectrum at wavenum- 
bers k < 0.052/i Mpc~^ is partly due to the effects of the 
uncertainty in the mean number density of optical galaxies 
(see Tadros & Efstathiou 1996 for a discussion of this ef- 
fect). However, independent evidence for the presence of a 
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Figure 1. The power spectrum of the distribution of galaxies and clusters of galaxies. Filled circles, open circles and crosses show the 
power spectrum of the galaxy distribution in the APM, Stromlo-APM and SSRS2-|-CfA2 surveys, respectively. Filled triangles and open 
squares show the power spectrum of the distribution of Abell clusters as determined by Einasto et al. (1997a) and Retzlaff at al. (1998), 
respectively. Open triangles represent the power spectrum of APM clusters. For comparison, we show the linear power spectrum of 
density fluctuations in the flat CDM models with Uq = 0.3 and h = 0.7 (solid line) and h = 0.6 (dashed line). The dotted line shows the 
power spectrum in the CDM model with Q = 1 and h = 0.5. The CDM models are COBE-normalized. 



preferred scale in the Universe at about 120h~^ Mpc comes 
from an analysis of the distribution of galaxy clusters. Fig- 
ure 1 shows the power spectrum of the distribution of the 
Abell-ACO clusters as determined by Einasto et al. (1997a) 
and Retzlaff et al. (1998), and the power spectrum of the 
APM clusters as measured by Tadros, Efstathiou & Dalton 
(1998). The power spectrum of the distribution of the Abell- 
ACO clusters has a well-defined peaJi at the same wavenum- 
ber, ko = 0.052/1 Mpc~^, as the power spectrum of galaxies 
in the Stromlo-APM survey. For wavenumbers k > ko, the 
shape of the clusters' power spectrum is similar to the shape 
of the power spectrum for galaxies in the Stromlo-APM sur- 
vey. This comparison suggests that the peak observed in the 
power spectrum of the Stromlo-APM redshift survey is a real 
feature in the distribution of galaxies (see Gramann (1998) 
for a more detailed discussion of the observed power spectra 
in different galaxy surveys). 

Cosmological models based on coUisionless dark mat- 
ter (e.q. cold dark matter (CDM)) and adiabatic fluctua- 
tions, when combined with power-law initial power spec- 
tra, predict smooth power spectra of density fluctuations at 
z ~ 10^. Figure 1 shows the power spectra of density fluc- 



tuations predicted in the flat CDM models with the density 
parameter fio = 0.3 and the normalized Hubble constant 
h — 0.6 and h — 0.7. For comparison, we show in Figure 1 
the power spectrum predicted in the CDM model with f2 = 1 
and h = 0.5. We have used the transfer function derived by 
Bardccn et al. (1986) and Sugiyama (1995), and the COBE 
normalization derived by Bunn and White (1997). The ob- 
served power spectra of galaxies and clusters of galaxies are 
not consistent with CDM-type models (see also e.g. Pear- 
cock 1997; Einasto et al. 1997a; Tadros, Efstathiou & Dalton 
1998). The baryonic acoustic oscillations in adiabatic mod- 
els may explain the observed excess only if currently favored 
determinations of cosmological parameters arc in substan- 
tial error (e.g. the density parameter Qo < 0.2h) (Eisenstein 
et al. 1998). One possible explanation for the presence of a 
peak in the power spectrum is an inflationary scenario with 
a scalar field whose potential has a localized feature around 
some value of the field (Starobinsky 1992; Lesgourgues, Po- 
larski & Starobinsky 1998). 

In this paper we study the properties and spatial distri- 
bution of galaxy clusters in two cosmological models which 
start from the observed power spectra of the distribution of 
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Figure 2. The linear and non-lineaj- power spectra of density fluctuations in the model (1) (solid lines) and in the model (2) (dot-dashed 
lines). The heavy curves show the results of the N-body simulations and the light curves the corresponding linear power spectra. The 
dashed line shows the linear power spectrum derived by Peacock (1997). Filled circles, open circles and crosses show the power spectrum 
of the galaxy distribution in the APM, Stromlo-APM and SSRS2+CfA2 surveys, respectively. 



galaxies. In the first model we assume that the initial linear 
power spectrum of density fluctuations in the universe at 
z ~ 10^ is of the form 



P{k) 



Piko){k/h>) 



ii k < ko; 
if > feo 



(1) 



where ko = 0.052/i Mpc"^ and P(fco) = 3.71 x IO^ctI/i"^ 
Mpc^. The as is the rms mass fluctuation on the 8h~^ Mpc 
scale. Gramann (1998) used the function (1) to recover the 
power spectrum in the Stromlo-APM redshift survey. In the 
second model we assume that the linear power spectrum 
contains a primordial feature at wavenumbers k ~ 0.1 — 0.2/i 
Mpc~^ (A ~ 30 — 60/i~^ Mpc) which correspond to the scale 
of superclusters, and 



P{k) = < 



P(fco)(fc/fco), if k < fc„; 

P(/co)(A;/A;o)"^ if feo < fe < fei 

P{ki), ifki <k < /C2; 

P{ki){k/ko)-^, iik>k2, 



(2) 



where ko = 0.052/i Mpc"\ fci = Q.lh Mpc~\ = 0.2ft 
Mpc'S P(fco) = 3.34 X 10*o-|/i-^ Mpc^ and P{ki) = 
4.07 X 10^<t|/i~^ Mpc^. We assume also that the initial den- 
sity fluctuation field in the Universe is a Gaussian field. In 



this case the power spectrum provides a complete statistical 
description of the field. 

Figure 2 demonstrates the linear and nonlinear power 
spectra of density fluctuations in the model (1) for as = 0.8 
and in the model (2) for as = 0.84. To examine the nonlin- 
ear evolution of density fluctuations we have used N-body 
simulations. The parameters of the simulations are given in 
Section 3. We see that at wavenumbers k > 0.2h~^ Mpc, the 
nonlinear power spectra in models (1) and (2) are very simi- 
lar. Figure 2 shows also the power spectra of the galaxy clus- 
tering in the APM, Stromlo-APM and SSRS2+CfA2 surveys 
(Baugh & Efstathiou 1993; Tadros & Efstathiou 1996; Costa 
et al. 1994) . The spatial distribution of galaxies in the mod- 
els (1) and (2) depends on the relation between galaxies and 
matter density, and further study is needed to study galax- 
ies in these models. Figure 2 demonstrates that if we assume 
a linear bias between galaxies and density, the inital power 
spectra given by equations (1) and (2) are in good agreement 
with the observed power spectra of galaxies. For comparison, 
we show in Figure 2 the linear power spectrum derived from 
the observed galaxy power spectra by Peacock (1997) (we 
have used his eq. [34] with parameters given in eq. [35]). The 
initial power spectrum derived by Peacock (1997) is similar 
to the function (1). 
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In this paper we examine the evolution of the mass func- 
tion, peculiar velocities, the power spectrum and the corre- 
lation function of gala^xy clusters in models (1) and (2) for 
different values of the density parameter fio and erg ■ To study 
the mass function of clusters of galaxies, wo use the Prcss- 
Schechter (1974) formalism. To investigate peculiar veloci- 
ties and the spatial distribution of galaxy clusters, we use 
N-body simulations. 

In the simulations we use the equations of motion for 
the model with 0=1. According to the equations of mo- 
tion, expressed in terms of the linear growing mode, the 
evolution of a pressureless fluid in an expanding universe 
is almost independent of the density parameter Qq and of 
the cosmological constant A (e.g. Gramann 1993; Nusser & 
Colberg 1998). Nusser & Colberg (1998) used high resolu- 
tion N-body simulations to investigate the effect of changing 
the cosmological background on the evolution of fluctuations 
and demonstrated that once the initial density fluctuation 
field is evolved to a given amplitude (e.g. as ~ 0.7) and 
smoothed on scales R > lh~^ Mpc, it is almost insensitive 
to the cosmological background. The smoothed nonlineax 
velocity field scales with the linear velocity growth factor, 
/(fio) ~ f^o '', .just as it docs in the linear theory. Therefore, 
if clusters represent the max;ima of the density field which 
is smoothed on scales ~ 1 — 2/i~^ Mpc, their spatial dis- 
tribution in real space is not sensitive to f2o and A. 

Observations provide the distribution of clusters in the 
rcdsliift space, which is distorted due to peculiar velocities of 
clusters. In order to study peculiar velocities of galgixy clus- 
ters and their distribution in the redshift space in the models 
with different flo, we determine the velocities of clusters in 
the simulations with = 1 and assume that peculiar ve- 
locities of galaxy clusters, as the whole velocity field, are 
proportional to the growth factor /(Oo). The linear velocity 
growth factor depends very weakly on the cosmological con- 
stant (e.g. Lahav et al. 1991) and in this paper we neglect 
this dependence. 

This paper is organized as follows. In Section 2 we study 
the mass function of clusters of galaxies in our models and 
compare the results with observations. In Section 3 we ex- 
amine peculiar velocities of galaxy clusters. In Section 4 we 
investigate the redshift-spacc power spectrum of clusters and 
in Section 5 we study the correlation function of clusters. 
Section 6 summarizes the main results. 

A Hubble constant of Ho = 100ft km s~^Mpc~^ is used 
throughout this paper. 



between M and M + dM is given by 



2 THE MASS FUNCTION OF CLUSTERS OF 

GALAXIES 

To investigate the mass function of clusters we use the Press- 
Sclicchter (1974, PS) approximation. The PS mass function 
has been compared with N-body simulations (Efstathiou et 
al. 1988; White, Efstathiou & Prenk 1993; Laccy & Cole 
1994; Eke, Colo & Fronk 1996; Borgani et al. 1997a) and 
has been shown to provide an accurate description of the 
abundance of virialized halos of cluster size. In the PS ap- 
proximation the number density of clusters with the mass 



n{M)dM = - 



2 p6 St da{M) 
TrMa^iM) dM 



5? 



(3) 

Here pb is the moan background density and St is the linear 
theory overdonsity for a uniform spherical fiuctuation which 
is now collapsing; St = 1.686 for O = 1, with a weak depen- 
dence on Oo for flat and open models (e.g. Eke et al. 1996). 
The function a{M) is the rms linear density fiuctuation at 
the mass scale M. We will use the top-hat window function. 
For the top-hat window, the mass M is related to the win- 
dow radius i? as M = AwpbR^/S. In this case, the number 
density of clusters of mass larger than M can be expressed 



nci(>M)= / n{M')dM' 



(27r)3/2 



St da{r) 
a^{r) dr 



exp 



2a2(-r) 



dr 



(4) 



Figure 3 shows the cluster mass function for the power 
spectra (1) and (2). We investigated the cluster masses 
within a 1.5h~^ Mpc radius sphere around the cluster cen- 
ter. This mass Mi. 5, is related to the window radius R as 



R = 8A3n^- 



Mi 



6.99 X W^^Qoh-^Me 



(/i"'Mpc). 



(5) 

Here the parameter a describes the cluster mass profile, 
M(r) ^ r" , at radii r ~ 1.5h~^ Mpc. Numerical simula- 
tions and observations of clusters indicate that the param- 
eter a ~ 0.6 — 0.7 for most of clusters (Navarro, Frenk & 
White 1995; Carlberg, Yee & EUingson 1997). In this paper 
we use a value a = 0.65. 

The number density of massive clusters is a sensitive 
function of Oq and ag- Figure 3a shows the mass function 
for the open models with Oo = 0.3, as = 0.8 and as = 1.0. In 
Figure 3b we present the mass function for the open mod- 
els with Oo = 0.4, as = 0.7 and as = 0.9. The cluster 
abundances in models (1) and (2), when compared at the 
same values of Oo and as, are very similar for smaller clus- 
ter masses. Mi. 5 < 4 x W^*h~^MQ. For larger masses the 
mass function in the model (2) is steeper than in the model 
(1). We investigated also the cluster abundances for the fiat 
models with slightly larger values of St derived by Eke et al. 
(1996), and found that for a given Oq and Mi. 5, the cluster 
abundance for the fiat model is ~ 10% smaller than for the 
open model. 

Figure 3 shows also the mass function of clusters of 
galaxies derived by Bahcall and Cen (1993, BC) and by Gi- 
rardi et al. (1998, G98). BC used both optical and X-ray 
observed properties of clusters to determine the mass func- 
tion of clusters. The function was extended towards the faint 
end using small groups of galaxies. G98 determined the mass 
function of clusters by using virial mass estimates for 152 
nearby Abell-ACO clusters including the new ENACS data 
(Katgcrt et al. 1998). The mass function derived by G98 is 
somewhat larger than the mass function derived by BC, the 
difference being larger at larger masses (see Figure 3). We 
find that the models with Oq = 0.3, as = 0.8 and fio = 0.4, 
as = 0.7 provide good match to the mass function derived by 
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Figure 3. The cluster mass function in the model (1) (solid lines) and in the model (2) (dot-dashed lines), (a) The density parameter 
Qq = 0.3. The mass function is shown for erg = 0.8 and for = 1.0. (b) The density parameter Qq = 0.4. The mass function is shown 
for erg = 0.7 and for erg = 0.9. Open circles and squares show the mass function of galaxy clusters derived by Bahcall and Cen (1993) 
and by Girardi et al. (1998), respectively. The open triangle describes the result obtained by White, Efstathiou & Prenk (1993). 



BC, whereas models with Oo = 0.3, ag = 1.0 and Qo = 0.4, 
as = 0.9 are in good agreement with the data derived by 
G98. 

Lot us consider the amplitude of the mass function of 
galaxy clusters at Mi. 5 = 4xl0^'^h~'^ Mq. For this mass, the 
cluster abundances derived by BC and G98 are ri(> M) = 
(2.0±l.l)xl0-^/i^ Mpc"^ andn(> M) = (6.3±1.2)10-'^/i^ 
Mpc~^, respectively. By analysing X-ray properties of clus- 
ters. White, Efstathiou & Frenk (1993) found that the num- 
ber density of clusters with mass Mi. 5 w 4.2 x 10^*h~^MQ is 
n(> M) = 4 X W~^h^ Mpc~^. Figure 4 shows the hmits for 



asQ,Q , assuming that the mass function of galaxy clusters 
at Ml. 5 = 4 X lO"/i"^M0 is in the range (2 - 6.5) x lO'^h^ 
Mpc~^. These limits are similar in both our models. For 
fio = 0.3 and Qo = 0.4, we find that erg = 0.90 ± 0.12 and 
0-8 = 0.80 ± 0.09, respectively. For O = 1, cts = 0.56 ± 0.05. 
These limits for as are very similar to the limits derived by 
Eke et al. (1996) for the CDM models by analyzing X-ray 
temperatures of clusters. 
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3 PECULIAR VELOCITIES OF CLUSTERS OF 
GALAXIES 

To investigate peculiar velocities of galaxy clusters and their 
spatial distribution we used N-body simulations. The sim- 
ulations examined in this paper were created using the 
particle-mesh code described by Gramann (1988). We in- 
vestigated the evolution of 256'^ particles on a 256^ grid, the 
comoving box size was L = 384/1-1 Mpc. Clusters were de- 
termined in the simulations at the moment when as = 0.8 
and (78 = 0.84 in the models (1) and (2), respectively. 

Clusters were selected in the simulations as maxima of 



the density field that was determined on a 256^ grid using 
the CIC-scheme. To determine peculiar velocities of clus- 
ters, we determined the peculiar velocity field on a 256^ 
grid using the CIC-scheme and found the peculiar veloci- 
ties at the grid points were the clusters had been identified. 
The clusters were then ranked according to their density 
and we selected Nd = {L/dd)* highest ranked clusters to 
produce cluster catalogs with a mean intercluster sepaxa^ 
tion del = 10 — lOO/i-i Mpc. For comparison, the num- 
ber density of the observed APM clusters and Abell clus- 
ters is rid ~ 3.4 X IQ-^/i^ Mpc"^ {dd ~ Sl/i'^Mpc) and 
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Uci ~ 2.5 X lO'^h^ Mpc~^ {del ~ 34/i"^ Mpc), respectively 
(Dalton et al. 1994, Einasto et al. 1997b, Retzlaff et al. 1998). 

It is difficult to follow the evolution of rich clusters by 
using N-body simulations, as it requires simulations with a 
very large dynamical range to identify correctly a sufficient 
number of clusters. The grid size in our N-body simulations 
is Rg = 1.5h~^ Mpc and the cluster centers in the simular- 
tions are at least 3h~^ Mpc apart. This could cause clusters 
to merge prematurely, and to have some effect on their prop- 
erties. However, in order to increase the resolution of the 
simulations we must increase the number of test particles 
and grid points, or have to follow the evolution of clusters 
in a smaller box. While the first possibility is technically 
difficult, in the latter case the rmmbcr of rich clusters be- 
comes too small to get statistically reliable results. Taking 
into account the requirements on the number of clusters and 
on the resolution together with the fact of fixed computer 
resources we decided to use a box size L = 384/i~^ Mpc and 
a grid size Rg = 1.5h~^ Mpc. 

To determine the rms peculiar velocities of clusters, we 
used the equation 



2 

Vcl 



2 I 2 



1 V ^ 2 2 



(6) 



where the parameter describes the dispersion of the clus- 
ter velocities, Vi, derived from the simulation and the pa- 
rameter u_L is given by 



1 1 



fe2 



(7) 



The last term in this expression is a discrete sum over the 
linear modes in the simulation, ks, with the power per mode, 
Ps(k), as actually used in the simulation. We found that for 
clusters with a mean separation dd > 30fe~^ Mpc, v\/v1 ~ 
4.9% and vl/vl ^ 3.7%, in models (1) and (2), respectively. 

Figure 5 shows the rms peculiar velocities of clusters, 
Vcl, in models (1) and (2), for the same values of f2o and 
as as in Figure 3. The rms peculiar velocities of clusters in 
the model (2) are ~ 12% smaller than in the model (1), 
when compared at the same values of f2o and ug. Figure 5 
shows the rms peculiar velocities for different values of the 
mean cluster separation, dd- We see that the rms peculiar 
velocities of clusters increase with cluster richness. In both 
models studied, the rms peculiar velocity of very massive 
clusters with an intercluster separation dd = 80/i~^ Mpc is 
^ 25% higher than the rms velocity of the clusters with a 
separation dd = 2Qh~^ Mpc. 

We note that this effect is sensitive to the agorithms 
used to define clusters and to determine their peculiar ve- 
locities. When we determined the clusters using the standard 
friends-of-friends algorithm (FOF) and defined the peculiar 
velocity of each cluster to be the mean peculiar velocity of 
all the particles within the cluster, we found that the rms 
peculiar velocities of clusters decrease with cluster richness. 
This is because with the FOF algorithm we identify almost 
the same objects as with the algorithm used throughout this 
paper, but the sizes of objects are different poor clusters 
arc smaller and rich clusters arc larger. The peculiar velocity 
of the cluster, as the mass of the cluster, depends on the size 
of the cluster. The larger the smoothing region, the smaller 
the velocity. 



In this paper we have determined the clusters as the 
maxima of the density field smoothed on the scale R ~ 
1.5/i~^ Mpc and have defined their peculiar velocities us- 
ing the same smoothing scale as for the density field. In this 
way, the sizes of all clusters are the same and do not depend 
on the richness of the cluster. Our study shows that in this 
case the rms peculiar velocities of clusters increase with clus- 
ter richness. In other words, the rms peculiar velocities of 
peaks increase with the height of the peaks. We found that 
this result is not very sensitive to the power spectrum of 
density fluctuations. We tested also two other models with 
different power spectra (n = — 1 and a sCDM model) and 
found that the rms peculiar velocity increases with cluster 
richness similary in all models studied. 

Our result is consistent with the result found by Colberg 
et al. (1998) for superclusters. Richer clusters axe clustered 
more strongly and, therefore, in superclusters the rms mass 
of clusters is larger than for isolated clusters. Colberg et al. 
(1998) studied the peculiar velocities of clusters and found 
that the rms velocities for clusters which are members of 
superclusters are about 20% to 30% larger than those for 
isolated clusters. Therefore, more massive clusters in super- 
clusters move faster than less massive clusters outside su- 
perclusters. 

To examine the rms peculiar velocities of clusters for 
different moments between as = 0.7 to erg = 1.0, we started 
from the cluster velocities at the moment as = 0.8 and 

as = 0.84, in models (1) and (2), respectively, and used the 
linear scaling, Vd ^ as- We tested this scaling in the model 
(1), by comparing the cluster velocities in N-body simula- 
tions at the moments when as = 0.5 and as = 0.8. The re- 
sults are presented in Figure 5b. For Oo = 0.4 and as ~ 0.7, 
we show the rms cluster peculiar velocities determined by 
using N-body simulations of clusters when as = 0.5 (dashed 
line) and as = 0.8 (solid line). We see that the differences 
for clusters with dd < 80h~^ Mpc are very small. The dif- 
ference is somewhat larger for massive clusters with a mean 
separation dd > 80h~^ Mpc (~ 5%). Therefore, during the 
evolution between as = 0.5 and as = 0.8, cluster velocities 
evolve almost as expected by the linear approximation. 

The observed rms peculiar velocity of galaxy clusters 
was investigated in several recent papers (e.g. Bahcall, Gra- 
mann & Gen 1994, Bahcall and Oh 1996, Borgani et al. 
1997T), Watkins 1997). In this paper we use the results ob- 
tained by Watkins (1997). He developed a likelihood method 
for estimating the rms peculiar velocity of clusters from line- 
of-sight velocity measurements and their associated errors. 
This method was applied to two observed samples of clus- 
ter peculiar velocities: a sample known as the SCI sample 
(Giovanelli et al. 1997) and a subsamplc of the Mark III 
catalog (WiUick et al. 1997). Watkins (1997) found the rms 
one-dimensional cluster peculiar velocity of 265lyg^ km s~^, 
which corresponds to the three-dimensional rms velocity of 
459+}f,(5 km s-\ 

Figure 4 shows the limits for asfio'^ in difi'erent mod- 
els, assuming that the observed cluster sample studied by 
Watkins (1997) corresponds to the model clusters with a 
mean cluster separation dd ~ 30h~^ Mpc {rid = 3.70 x 
lO'-'h^ Mpc"-'*). In the model (1), the rms peculiar velocity 
of clusters with a separation dd = 30h~^ Mpc is 4591 J^gg 
km s-\ when as = {0.33t^ ll)f-^{no). For fio = 0.3 and 
fio = 0.4, we obtain as = O-GStoAl and as = 0.57lo i6, 
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Figure 5. The rms peculiar velocities of clusters in the model (1) (solid lines) and in the model (2) (dot-dashed lines). The rms velocities 
are shown for different values of the mean cluster separation, d^i- (a) The density parameter f2o = 0.3. The rms velocities are shown 
for (Tg = 0.8 and for erg = 1.0. (b) Tlie density parameter = 0.4. The rms velocities are shown for trg = 0.7 and for <Tg = 0.9. For 
(Tg = 0.7, we show the rms cluster velocities determined by simulations when eg = 0.5 (dashed line) and (Tg = 0.8 (solid line). The filled 
circle shows the observed rms peculiar velocity of galaxy clusters. 



respectively. For the power spectrum (2), we found that 
(Tg = (0.37lJ5:i^)/"^(t7(,). Therefore, for = 0.3 and 
Qo = 0.4, (Tg = 0.76l!!:fi and as = 0.64^',;;^?, respectively. 

Now we can compare the observational constraints ob- 
tained by studying the mass function and peculiar veloci- 
ties of clusters of galaxies. In the model (1) for fio = 0.4, 
the mass function and the peculiar velocities of clusters are 
consistent with the observed data, for the small window of 
(Tg = 0.71 — 0.78. For Qq > 0.6, the observed mass func- 
tion and the peculiar velocities of clusters are not consistent 
with each other. Either the observed mass function of clus- 



ters is overestimated or the peculiar velocities of clusters arc 
underestimated. Therefore, in the model (1), the mass func- 
tion and the peculiar velocities of clusters are consistent with 
observations only if Qq < 0.6. In the second model, the per- 
mitted window in the {flo, as) plane is larger. For Oo — 0.4, 
the mass function and the peculiar velocities arc consistent 
with the observed data if ag = 0.72 — 0.88. For example, 
for = 0.4 and as = 0.75, the number density of clusters 
with mass Mi. 5 = 4 x IO^/i-^Mq is n(> M) = 2.7 x IQ-^/i^ 
Mpc^"* , and the rms peculiar velocity of clusters with a mean 
separation dd = 30h~^ Mpc is 530 km s~^. 
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Figure 6. The rms peculicir velocity of peaJfs, ap{R), for the power spectra (1) (heavy solid line) and (2) (heavy dot-dashed line) for 
Ho = 0.3 and erg = 0.9. The light curves show the corresponding rms peculiar velocity crv{R) for the same models. The small and large 
squares demonstrate the rms peculiar velocity of clusters in the model (1) with mean separations dd = 10h~^ Mpc and dd = 30h~^ 
Mpc, respectively. The small and large circles show the corresponding rms peculiar velocities in the model (2). 



We investigated also the linear theory predictions for 
pecuUar velocities of peaks in the Gaussian field. The linear 
rms velocity fluctuation on a given scale R can be expressed 
as 

a4R) = f/o/(no)a-i(i?), (8) 
where aj is defined for any integer j by 

a| = ^ / P{kW''{kR)k^'+^dk. (9) 

Bardeen et al. (1986) showed that the rms peculiar velocity 
at peaks of the smoothed density field differs systematically 
from av{R), and can be expressed as 

cTp{R) = a4R)^yi-cTt/afcTl,. (10) 

In this approximation, the rms velocities of the peaks do not 
depend on the height of the peaks. 

Figure 6 shows the rms peculiar velocities of peaks, 
ap{R), for the power spectra (1) and (2) for fio = 0.3 and 
as = 0.9. We have used the top-hat window function. For 
comparison, wc show also the rms peculiar velocity av{R) 
for the same models. For the cluster radius R = l.5h~^ Mpc, 
Up is lower than cr„ about ~ 2% and ~ 3.5% for models (1) 
and (2), respectively. On larger scales, the difference between 



ffp and ay increases. For comparison, we show in Figure 6 
the rms peculiar velocity of clusters with mean separations 
del = Wh^^ Mpc and dd = 30h^^ Mpc, derived by using N- 
body simulations. The rms peculiar velocity of clusters with 
a separation dd = Wh~^ Mpc is similar to the linear the- 
ory expectations at the scale R ~ 1.5/i~^ Mpc. (It is slightly 
smaller (~ 2%) than ap at the radius R = 1.5h~^ Mpc. This 
small difi'erenco is probably caused by smoothing inherent 
to particle-mesh method). 

The peculiar velocity of rich clusters is higher than that 
predicted by the linear approximaton (10). In the model 
(1), the peculiar velocity of clusters with a mean separation 
dd ~ 30h~^ Mpc is ~ 19.0% higher than ap at the radius 
R = l.5h~^ Mpc. In the model (2), the peculiar velocity 
of these clusters is ~ 16.8% higher than that predicted by 
the linear theory. These results are in good agreement with 
the results obtained by Colborg et al. (1998). They com- 
pared the peculiar velocities of galaxy clusters with ap{R) 
at larger radii, 7? ~ 8h'^ Mpc, and for this reason found a 
larger difference (~ 40%) between the rms peculiar velocity 
of galaxy clusters and ap. 

In Figure 5, we showed that during the evolution be- 
tween as = 0.5 and as = 0.8, the rms peculiar velocities 
of clusters with a mean separation dd < 80h~^ Mpc evolve 
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Figure 7. Panel (a) shows the distribution of clusters in a 96/i~^ Mpc thick slice in the model (1). Panel (b) shows the distribution of 
clusters in the same slice in the model (2). The distribution is shown for the clusters with a mean separation d^i = 34/i~^ Mpc. 



as expected by the linear theory, Vd ~ era. The difference 
between the linear approximation (10) and the peculiar ve- 
locities of the del ~ 30/i~^ Mpc clusters must therefore arise 
when as < 0.5. Further study is needed to investigate the 
nonlinear evolution of peculiar velocities of clusters in more 
detail. 



4 THE POWER SPECTRUM OF CLUSTERS 
OF GALAXIES 

Figure 7 shows the spatial distribution of clusters with a 
mean separation dd — 34/i^^ Mpc in our models. Figure 7a 
shows the distribution of clusters for the initial power spec- 
tra (1) and Figure 7b for the initial spectra with a primordial 
feature at wavelengths A ~ 30 — 60h~^ Mpc (equation [2]). 
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Figure 8. The rcdshift-spacc power spectrum of clusters in the model (1) for erg = 0.8 (solid lines) and in the model (2) for erg = 0.84 
(dot-dashed lines). The heavy curves show the power spectra of clusters and the light curves the corresponding linear power spectra 
of matter fluctuations. To transform the clusters to the redshift space we have assumed that Hq = 0.3. The dotted line demonstrates 
the redshift-space power spectrum of clusters in the model (1) for U = 1. The dashed line represents the power spectrum of clusters 
in the model (1) for erg = 0.5. (a) The power spectrum of model clusters and APM clusters (open triangles) with a mean separation 
d^l = 31h~^ Mpc. (b) The power spectrum of model clusters and Abell clusters with a separation dcj, = 34ft.~^ Mpc. Filled triangles and 
open squares show the power spectrum of Abell clusters as determined by Einasto et al. (1997a) and RetzlaiT at al. (1998), respectively. 



The phases of the initial density fluctuations were chosen 
to be same, and therefore, we see directly the influence of 
the initial power spectrum. At a first look, the distribution 
of clusters in different models is rather similar. However, 
we see that the shape of superclusters in different models 
is sUghtly different. In the model (1), the clusters inside su- 
perclusters axe more concentrated, while in the model (2), 
superclusters are larger and the clusters within superclusters 



are more disperse (compare, e.g., the large superclusters in 
the upper-left corner). The bump in the initial power spec- 
trum at A ~ 30 — 60h~^ Mpc influences the distribution of 
clusters inside superclusters. 

Figure 8 demonstrates the redshift-space power spec- 
trum of clusters, Pci{k). The clusters were determined in the 
simulation at the moment when as = 0.8 and as = 0.84, in 
the models (1) and (2), respectively. To calculate the power 
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spectrum wc transformed the cluster positions to the red- 
shift space, determined the density field on a 128^ grid us- 
ing the CIC scheme and calculated its Fourier components, 
subtracting the shot noise term. 

To transform the clusters from the real space to the red- 
shift space, we determined the peculiar velocities of clusters 
in the simulations with = 1 and used the linear scaling, 
Vci ~ /(f^o), for fio = 0.3. In this case the peculiar veloci- 
ties of clusters are consistent with observations for <T8 ~ 0.8 
(see Figure 5). However, Pci{k) is not very sensitive to this 
assumption. In Figure 8 we show also the power spectrum 
of clusters in the redshift space for the model (1) for f2 = 1, 
where the peculiar velocities of clusters are severly overesti- 
mated, Vci ~ 1100 km s^^. In this model Pci{k) at wavenum- 
bers k < O.l/i"^ Mpc is - 10% higher than for fio = 0.3. 

We investigated also the cluster power spectrum in the 
model (1) for as = 0.5 (here we used n = 1). We found 
that during the evolution between erg = 0.5 and ug = 0.8, 
the power spectrum of clusters with a mean separation 
del < 30h~^ Mpc is almost unchanged. For richer clusters, 
Pci{k) somewhat decreases (~ 14% and ~ 29% for the clus- 
ters with del = Slh~^ Mpc and dd = 34/i~^ Mpc, respec- 
tively.) This effect is probably caused by merging of very 
rich clusters. Further study (e.g. numerical simulations with 
higher dynamical range) is needed to determine whether this 
is a real effect for the model (1), or a numerical effect due 
to the limited dynamical range of the N-body simulations. 

Let us compare the power spectra of clusters in our 
models with the observed power spectra of the APM and 
Abell clusters. Figure 8a shows the power spectrum of model 
clusters with a moan separation dd — 31h~^ Mpc. For com- 
parison, we show the power spectrum of the observed APM 
clusters determined by Tadros, Efstathiou & Dalton (1998). 
They analyzed the redshift survey of 364 clusters described 
by Dalton et al. (1994). The mean intercluster separation of 
the APM clusters is d^i ~ 31/r^Mpc (Dalton ct al. 1994). 
Figure 8a shows that the power spectrum of the APM clus- 
ters is in good agreement with the power spectrum of clus- 
ters predicted in the model (2). In the model (1), the power 
spectrum of clusters is higher than observed (factor of ~ 1.8 
at fe ~ 0.07 - 0.08/1 Mpc^i). 

Figure 8b demonstrates the power spectrum of model 
clusters and the Abell clusters with a mean separation 
del = 34/i~^ Mpc. We show the power spectra of the Abell 
clusters determined by Einasto et al. (1997a) and Retzlaff 
ct al. (1998). Einasto ct al. (1997a) determined the power 
spectrum of the Abell clusters from the correlation function 
of clusters, while Retzlaff et al. (1998) estimated the power 
spectrum directly from the Fourier space. Wc sec that the 
power spectra estimated by different methods arc consistent 
with each other. However, the error bars measured by Ret- 
zlaff et al. (1998) are much larger than calculated by Einasto 
et al. (1997a) and probably underestimated in the latter case 
(see Retzlaff et al. 1998, Einasto et al. 1997a for details). 
Figure 8b shows that the power spectrum predicted in the 
model (1) for erg = 0.8 is consistent with the observed power 
spectrum of the Abell clusters. For erg = 0.5, the amplitude 
of density fluctuations at wavenumbers k > O.lh Mpc~^ is 
higher than observed. The power spectrum of clusters in the 
model (2) is consistent with the observed power spectrum 
measured by Retzlaff et al. (1998) within the uncertainties. 

We investigated also the relation between the power 



spectrum of clusters and the power spectrum of matter 
fluctuations. During the evolution the power spectrum of 
clusters, Pci{k), is almost unchanged, while the power spec- 
trum of matter fluctuations evolves as P{k) ~ erf in the 
linear regime. In this case, the bias parameter b^iik) = 
Pci{k)/P{k) ^ and therefore, in order to compare differ- 
ent models at different moments, it is reasonable to express 
the parameter bd in terms of erg. We examined the bias pa- 
rameter bd at the wavenumber interval k ~ 0.06 — O.OSh 
Mpc^^. In the model (1) for ag — 0.8, wc found that 
bd = 2.30/(Tg and bd = 2.40/CTg for the clusters with 
dd = 31h~^ Mpc and dd = 34ft~^ Mpc, respectively. (For 
(Tg = 0.5, as the cluster power spectrum is somewhat larger, 
the bd = 2.45/(78 and bd ~ 2.70/(Tg, respectively). In the 
model (2) for a — 0.84, we obtain that bd = 2.25/crg and 
bd = 2.50/ as, respectively. In this model, the bias param- 
eter is similar to the model (1). The relation between the 
distribution of matter density and of clusters of different 
richness is studied in more detail in the paper by Gramann 
& Suhhonenko (1998). 



5 THE CORRELATION FUNCTION OF 
CLUSTERS OF GALAXIES 

In this Section we examine the spatial two-point correlation 
function of clusters, Cci(^)) our models. To determine the 
correlation function, we first tranformed cluster positions 
to the redshift space and then determined the correlation 
function of clusters by pair counting. 

Figure 9 shows the redshift-spacc correlation function of 
clusters in the model (1) for erg = 0.8 and in the model (2) 
for erg = 0.84. The peculiar velocities were calculated for the 
density parameter Qg = 0.3. We investigated also the f2 = 1 
model for velocities and found that in this case ^d{r) at 
radii r = 10-20/r^Mpc is - 16% higher than for Qo = 0.3. 
Figure 9 shows also the correlation function in the model 

(1) for as ~ 0.5. During the evolution in the model (1), the 
correlation function of clusters, as the power spectrum of 
clusters, somewhat decreases (~ 22% and ~ 38%, at r ~ 
15/i"^ Mpc, for the dd = 31h'^ Mpc and dd = 34/i"^ 
Mpc clusters, respectively). This effect is probably caused 
by merging of rich clusters. 

Now we can compare the correlation function of clusters 
predicted in the models with the observed correlation func- 
tion of the APM and Abell clusters. Figure 9a demonstrates 
the correlation function of model clusters and of the APM 
clusters with a mean separation dd = Slh~^ Mpc. The cor- 
relation function of the APM clusters has been determined 
by Dalton et al. (1994). On small scales the correlation func- 
tion in the model (2) is in good agreement with the correla- 
tion function of the APM clusters. The correlation function 
of the APM clusters is equal to unity at a pair separation 
ro = 14.3 ± 1.75h~^ Mpc (Dalton et al. 1994). In the model 

(2) we find that ro = 14 ± Ih'^ Mpc. 

On larger scales, ^d{r) in the model (2) is larger than 
that measured by Dalton et al. (1994) for the APM clusters. 
On the other hand, the power spectrum of clusters on large 
spatial scales in the model (2) is in good agreement with the 
power spectrum of the APM clusters determined by Tadros, 
Efstathiou & Dalton (1998) (Figure 8a). This comparision 
suggests that the correlation function of clusters determined 
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Figure 9. The rcdshift-spacc correlation function of clusters in the model (1) for as = 0.8 (solid lines) and in the model (2) for 
erg = 0.84 (dot-dashed lines). The heavy curves show the correlation function of clusters and the light curves the corresponding linear 
correlation function of matter fluctuations. To transform the clusters to the redshift space we have assumed that = 0.3. The dotted 
line demonstrates the redshift-space correlation function of clusters in the model (1) for H = 1. The dashed line represents the correlation 
function of clusters in the model (1) for as = 0.5. (a) The correlation function of model clusters and APM clusters (open triangles) with 
a mean separation dd = 31h~^ Mpc. (b) The correlation function of model clusters and Abell clusters (filled triangles) with a separation 
del = 34/1-1 Mpc. 



by Dalton et al. (1994) may bo underestimated at large sep- 
arations r > 20h~^ Mpc. In the model (1), i,ci{r) is higher 
than observed on all scales. In this model ro = 19 ± lh~^ 
Mpc. 

Figure 9b shows the spatial correlation function of clus- 
ters with a separation dd = 34/i~^ Mpc. For comparison, 
we present the correlation function of the Abell clusters de- 
termined by Einasto et al. (1997c). On small scales the cor- 

© 0000 RAS, MNRAS 000, 000-000 



relation function of the Abell clusters is higher than that 
observed for the APM clusters. We find that ro = 20±3/i"^ 
Mpc for the Abell clusters. This effect is partly due to dif- 
ferences in the number densities of the APM and Abell clus- 
ters, liowcvcr, our models predict that the correlation func- 
tion of clusters with a mean separation dd = 31h~^ Mpc 
and dd = 34/i.~^ Mpc is very similar. (In the models (1) 
and (2), the correlation length of clusters with dd = 34/1" ^ 
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Mpc is ro = 20 ± lh~^ Mpc and ro = 15 ± Ih'^ Mpc, re- 
spectively.) The correlation function of the Aboll clusters 
at separations r < 25h~^ Mpc is probably overestimated 
due to the projection and selection biases known to affect 
clustering in the Abell cluster catalogues (e.g. Sutherland 
1988). On larger separations these effects are not so impor- 
tant. The correlation function of clusters in the models (1) 
and (2) is consistent with the observed correlation function 
of the Abell clusters on separations r > 25h~^ Mpc. We 
examined also the radius at which the cluster correlation 
function ^ci{ri) = 0. We found that n = 54 ± 3h~^ Mpc 
and ri = 60 ± 3h~^ Mpc, in the models (1) and (2), respec- 
tively. For the Abell clusters, the parameter n = 50±10/i~^ 
Mpc. 

Thus, the model (2) fits the correlation function of the 
APM clusters on small scales and the correlation function of 
Abell clusters on large scales. This model predicts that there 
is a bump in the correlation function of clusters at separa- 
tions r ~ 20 — 35/i^^Mpc. Available data are insufficient to 
confirm or to rule out this interesting possibility. Accurate 
measurements of the correlation function of clusters at these 
distances can serve as a discriminating test for this model. 



6 CONCLUSIONS 

In this paper, we have examined the properties of clusters of 
galaxies in two cosmological models. In the first model, the 
initial power spectrum was chosen in the form P{k) oc at 
the wavelengths A < 120/i~^ Mpc (equation 1). In the second 
model, we assumed that the initial power spectrum contains 
a primordial feature at the wavelengths A ~ 30 — 60/i^^ Mpc 
(equation 2). The density fluctuations at these wavelengths 
influence the distribution of clusters inside superclusters. In 
the model (2), superclusters are larger and the clusters inside 
superclusters arc not so concentrated than in the model (1) 
(see Figure 7). We investigated the mass function, peculiar 
velocities, the power spectrum and the correlation function 
of clusters in both models for different values of f2o and as- 
Below, we briefly summarize the results obtained. 

(1) The mass function of clusters of galaxies in models 

(1) and (2), when compared at a same values of Oq and as, 
is very similar for smaller masses M < 4 x 10^'^h~^ Mq. For 
larger masses the mass function in the model (2) is steeper 
than in the model (1). For fio = 0.3 and Qo = 0.4, the 
mass function of clusters in both our models is consistent 
with observations, if aa = 0.90 ± 0.12 and as = 0.80 ± 0.09, 
respectively. 

(2) The rms peculiar velocities of clusters in the model 

(2) are ~ 12% smaller than in the model (1), when com- 
pared at the same values of ilo and as- In the model (1), 
the rms peculiar velocity of clusters is consistent with ob- 
servations if (78 = (0.33l" (',fj)n,7'''®. In this model, the mass 
function and the peculiar velocities of clusters are consistent 
with observations only if flo < 0.6. For ilo = 0.4, the mass 
function and the peculiar velocities are consistent with the 
observed data if as = 0.71 — 0.78. In the model (2), the rms 
peculiar velocity of clusters is consistent with observations 
if CTg = (0.37lo' j"(',)0(7'' '' and the permitted region in the 
{Qo,as) plane is larger. For f2o = 0.4, the mass function and 
the peculiar velocities are consistent with the observed data 
if 0-8 = 0.72 - 0.88. 



(3) The redshift-space power spectrum of clusters in the 
model (2) is in good agreement with the observed power 
spectrum of the APM clusters. The power spectrum of 
clusters in this model is also consistent with the observed 
power spectrum of the Abell clusters within uncertainties. 
In the model (1), the power spectrum of clusters is higher 
than observed for the APM clusters (factor of ~ 1.8 at 
fe ~ 0.07 - 0.08/1 Mpc"^). 

(4) The redshift-space correlation function of clusters 
in the model (2) is consistent with the correlation function 
of the APM clusters at small distances r < 25ft~^ Mpc. 
At larger separations the cluster correlation function in this 
model is consistent with the correlation function as derived 
for the Abell clusters. In the model (1), the correlation func- 
tion of clusters on small distances is higher than observed 
for the APM clusters. 

Therefore, in many aspects the power spectrum of den- 
sity fluctuations in the model (2) fits the observed data bet- 
ter than the simple power law model (1). The superclusters 
of galaxies in the Universe are probably more disperse as 
predicted in the model (2) and not so concentrated as pre- 
dicted in the model (1). Observed data are not sufficient 
to examine the power spectrum of density fiuctuations at 
wavelengths A ~ 30 - 120/i"^ Mpc in more detail, but our 
study suggests that probably at these wavelengths the initial 
power spectrum is not a featureless simple power law. 

We examined also the linear theory predictions for pe- 
culiar velocities of peaks in a Gaussian field and compared 
these with the peculiar velocities of clusters in N-body sim- 
ulations. We determined the clusters as the maxima of the 
density field smoothed on the scale R ~ 1.5h~^ Mpc and 
defined their peculiar velocities using the same smoothing 
scale as for the density field. In this way, the sizes of all 
clusters are the same and do not depend on the richness of 
the cluster. Our study shows that in this case the rms pecu- 
liar velocities of clusters increase with cluster richness. The 
rms peculiar velocity of small clusters is similar to the linear 
theory expectations, while the rms peculiar velocity of rich 
clusters is higher than that predicted in the linear theory. 
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